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The ground and some excited states of the Li atom in external uniform magnetic fields are 

calculated by means of our 2D mesh Hartree-Fock method for field strengths ranging from zero up 

CN ! to 2.35 ■ 10 8 T. With increasing field strength the ground state undergoes two transitions involving 

three different electronic configurations: for weak fields the ground state configuration arises from 

the field-free ls 2 2s configuration, for intermediate fields from the ls 2 2p_i configuration and in high 

fields the Is2p_i3d_2 electronic configuration is responsible for the properties of the atom. The 

transition field strengths are determined. Calculations on the ground state of the Li + ion allow us 

to describe the field-dependent ionization energy of the Li atom. Some general arguments on the 

ground states of multi-electron atoms in strong magnetic fields are provided. 
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I. INTRODUCTION 



The behaviour and properties of atoms in strong magnetic fields is a subject of increasing interest. On the o.h.s. this 
is motivated by the astrophysical discovery of strong fields on white dwarfs and neutron stars and on the o.h.s. 
the competition of the diamagnetic and Coulombic interaction causes a rich variety of complex properties which are of 
interest on their own. Investigations on the electronic structure in the presence of a magnetic field appear to be quite 
complicated due to the mixed geometry of this quantum problem (mixing of spherical and cylindrical symmetry). 
There are many works on the hydrogen atom (for a list of references see ) and several works on the He atom as 
well as He- like ions . Other atoms however have been investigated only in a very few cases [pT 13 lif]. 



For the hydrogen atom the impact of the mixed symmetry is particularly evident and at the same time pronounced 
in the intermediate field regime for which the magnetic and Coulomb forces are comparable. For different electronic 



degrees of excitation of the atom the intermediate regime is met for different absolute values of the field strength. 

•i-H ' 

' For the ground state the boundaries of this regime can be defined in a rough manner as the range 7 = 0.2 — 20 
(7 = B/Bq, B is the magnetic field strength, B — hc/eofi — 2.3505-10 5 T; atomic units will be used in the following). 
Oh With increasing degree of excitation the domain of the intermediate fields lowers correspondingly and becomes, as a 
rule, wider on a logarithmic scale of 7. Both early flql and more recent works |5|,|l6|-|l9| on the hydrogen atom have 
used different approaches for relatively weak fields (the Coulomb force prevails over the magnetic force) and for very 
strong fields where the Coulomb force can be considered as weak in comparison with the magnetic forces (adiabatic 
limit). In early works the Coulomb field was considered in this limit actually as perturbation for a free electron in 
a superstrong magnetic field. The motion of an electron parallel to the magnetic field is governed in the adiabatic 
approximation by a ID quasi-Coulomb potential with a parameter, dependent on the magnetic field strength. 
The detailed calculations of the hydrogen energy levels carried out by Rosner et al B also retain the separation of 
the magnetic field strength domains due to decomposing the electronic wave function in terms of either spherical 
(for weak fields) or cylindrical harmonics (for strong fields). A powerful method to obtain comprehensive results on 
low-lying energy levels in the intermediate regime in particular for the hydrogen atom is provided by mesh methods 
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For atoms with several electrons the problem of the mixed symmetries is even more intricate than for hydrogen 
because different electrons feel very different Coulomb forces, i.e. possess different single-particle energies, and the 
domain of the intermediate fields therefore appears to be the sum of the intermediate domains for the separate 
electrons. 

There exist several investigations on two-electron atoms in the literature f^- [f2] , [l^ , ^f| -p5| . The majority of them 
deals with the adiabatic limit in superstrong fields. Most of the early works are Hartree-Fock (HF) calculations for 
the strong field domain. There are also several variational calculations for the low-field domain (2^]2(||27j including 
calculations by Larsen |22| made at 7 < 2 for He atom and at 7 < 5 for H~. The latter calculations can be used for 
evaluations of the correlation energy in the low-field domain. HF calculations [^) are carried out analogously to the 
approach in ref. || with applying two different sets of basis functions to the high- and low-field domains. As a result 
of the complicated geometry in the intermediate regime this approach inherently suffers from very slow convergence 
properties with respect to the energy eigenvalues and yields therefore only a very low accuracy. Accurate calculations 
for arbitrary field strengths were carried out in refs. ||,[l(J by the 2D mesh HF method. Investigations on the ground 
state as well as a number of excited states of helium including the correlation energy have very recently been performed 
via a Quantum Monte Carlo approach Very recently benchmark results with a precision of 10~ 4 — 10~ 6 for the 
energy levels have been obtained for a large number of excited states with different symmetries using a configuration 
interaction approach with an anisotropic Gaussian basis set |2^] . 

For the lithium atom which is the subject of the present work there exists only one recent investigation by Jones et 
al [ |i"f| ] . It contains calculations for the ground state and a few low-lying states of the Li atom at weak and intermediate 
fields. Precise Hartree-Fock results for several states in weak fields and quite satisfactory results for the intermediate 
region are presented in this work. However their basis functions did not allow to perform calculations for stronger 
fields. An attempt to define the sequence of the electronic ground state configurations which are different for different 
regimes of the field strength has also been undertaken in this work. However a detailed qualitative analysis of the high- 
field ground state configuration was not carried out. As a result the high-field ground state electronic configuration 
and the transition point to this configuration from the intermediate one is still an open question. 

In the current work we apply a fully numerical 2D Hartree-Fock method to the problem of the Li atom in magnetic 
fields of arbitrary strength. This method enables us performing calculations for various states and with approximately 
equal precision for weak, intermediate and superstrong magnetic fields. Our main focus is the ground state of the Li 
atom and its ionization energies. To this end several electronic configurations of the Li atom and two configurations 
of the Li + ion are studied. 



II. FORMULATION OF THE PROBLEM AND METHOD OF SOLUTION 

We solve the electronic Schrodinger equation for the lithium atom in a magnetic field under the assumption of 
an infinitely heavy nucleus in the (unrestricted) Hartree-Fock approximation. The solution is established in the 
cylindrical coordinate system (p, <fi, z) with the z-axis oriented along the magnetic field. We prescribe to each electron 
a definite value of the magnetic quantum number m M . Each single-electron wave function vDy depends on the variables 
(j) and (p, z) 

& z) = (2n)- 1 ' 2 e-* m >+il> lt (z ) p) (1) 

where p = 1,2, 3 is the numbering of the electrons. The resulting partial differential equations for ip^z, p) and the 
formulae for the Coulomb and exchange potentials have been presented in ref. [jlO| . 

The one-particle equations for the wave functions p) are solved by means of the fully numerical mesh method 

described in refs. [f|[l0|]. The new feature which distinguishes the present calculations from those described in ref. 
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is the method of calculation of the Coulomb and exchange integrals. In the present work as well as in ref. 
we obtain these potentials as solutions of the corresponding Poisson equations. The problem of the boundary 
conditions for the Poisson equation as well as the problem of simultaneously solving Poissons equations on the same 
meshes with Schrodinger-like equations for the wave functions ipfi(z, p) have been discussed in ref. In the present 
approach these problems are solved by using special forms of non-uniform meshes. Solutions to the Poisson equation 
on separate meshes contain some errors dp associated with an inaccurate description of the potential far from the 
nucleus. However due to the special form of the function Sp(h) for these meshes (where h is a formal mesh step) the 
errors do not show up in the final results for the energy and other physical quantities, which we obtain by means of 
the Richardson extrapolation procedure (polynomial extrapolation to h — [|6|p9|] ) . An additional improvement with 
respect to the precision of our numerical calculations of the integrals is achieved by solving the Poisson equation not 
for the whole charge distribution but for the total distribution minus some properly chosen charge distribution with 
known analytical solution to the Poisson equation. Both of these approaches will be described in detail in a separate 
work. 

Our mesh approach is flexible enough to yield precise results for arbitrary field strengths. Some minor decrease of 
the precision appears in very strong magnetic fields. This phenomenon is due to a growing difference in the binding 
energies es M of single-electron wave functions belonging to the same electronic configuration 

(m„ + K,| + 2*,„ + 1)7/2 - e M (2) 

where e M is the single-electron energy and s zfl is the spin z-projection. This results in big differences with respect to the 
spatial extension of the density distribution for different electrons. This difference is important for the configurations 
ls 2 2s, Is2s2p_i and ls2po2p_i and is not important for Is2p„i3d_2 and ls 2 2p_i because all the single-electron 
energies for the latter states are of the same order of magnitude. The precision of our results depends, of course, on 
the number of mesh nodes and can be improved in calculations with denser meshes. The most dense meshes which 
we could use in the present calculations had 120 x 120 nodes. These meshes were used for the states ls 2 2s, Is2s2p_i 
and ls2po2p_i at fields 7 = 500 and 7 = 1000. For other states and weaker magnetic fields Richardson's sequences 
of meshes with maximal number 80 x 80 or 60 x 60 were sufficient. 

III. GROUND STATE ELECTRONIC CONFIGURATIONS 

We start this section with a qualitative consideration of the problem of the atomic multi-electron ground states 
in the limit of strong magnetic fields. It is clear that the state ls 2 2s of the lithium atom is the ground state only 
for relatively weak fields. The set of single-electron wave functions for constructing the HF ground state for the 
opposite case of extremely strong magnetic fields can be determined as follows. The nuclear attraction energies and 
HF potentials (which determine the motion along z axis) are then small compared to the interaction energies with the 
magnetic field (which determines the motion perpendicular to the magnetic field and is responsible for the Landau 
zonal structure of the spectrum). Thus, all the single-electron wave functions must correspond to the lowest Landau 
zones, i.e. < for all the electrons, and the system must be fully spin-polarized, i.e. s 2M = — \ (J.). For the 
Coulomb central field the single-electron levels form quasi ID Coulomb series with the binding energy Ep = -^1 for 
n z > and Ep — > 00 for n z — 0, where n z is the number of nodal surfaces of the wave function, which cross the z 
axis. These relations between single-electron energies and the geometry of single-electron wave functions along with 
analogous relations for the field-free atom provide the basis for the following considerations. 

It is evident, that the wave functions with n z = have to be choosen for the ground state at 7 — > 00. Thus, for 
7 — > 00 the ground state of the Li atom must be Is j 2p_i J, 3tf_2 j- This state was not considered in JIl[ but only 
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the Is I 2po I 2p_i I configuration was presented. Analogously, the very high-field ground state for the C atom 
considered in O] must be the state belonging to the configuration Is j 2p_i j 3cL 2 I 4/_3 J. 5g_4 J, 6/i_5 L 

The problem of the configuration of the ground state for the intermediate field region cannot be solved without doing 
cxplicite calculations. Calculations in ref. [fj"l[ were carried out for configurations with the maximal single-electron 
principal quantum number n < 2. Under this restriction calculations for the states ls 2 2s, ls 2 2p_i, Is J. 2s J, 2p_i j, 
and Is I 2po J, 2p_i j are sufficient to determine the set of intermediate ground states. Indeed, ls 2 2s is the zero-field 
ground state. ls 2 2p_i is the lowest excited state of the field free atom and (contrary to ls 2 2s) all the single-electron 
wave functions of this state must have infinite binding energies in the infinite strong magnetic field. Moreover, this 
state has the largest binding energy Eb 

3 

E B = ^(77V+|77V|+2s ZM + 1)7/2 -E (3) 

in the strong field limit due to the fact that es(ls) > es(2p_i) > cb(3c£_2) > ... in strong fields. (For 7 = 1000 one 
can obtain binding energies from table | as Eb(Is 2 2p-i) = 69.1569 and £ l s(ls2p_i3d_2) = 60.0589). The reader 
should note that the ls 2 2p_i configuration cannot represent the ground state in very strong fields since it is not fully 
spin polarized. The state Is J, 2s J, 2p_i J, is the lowest fully spin-polarized state with the single-electron principal 
quantum numbers < 2 in weak fields and, at last, the state Is J, 2p$ j 2p_i j which lies higher at 7 = must 
become lower than Is J, 2s J, 2p_i J, with increasing field strength. 

Our calculations include the high-field ground state Is J, 2p_i j 3(i_2 I which contains one electron with n = 3. In 
principle, also other configurations could be considered as possible ground states for intermediate field strength. Such 
configurations are ls 2 3s, ls 2 3p_i, ls 2 3d_2, Is2s3s, Is2s3p_i, Is2s3d_2, Is2p_i3s, and Is2p_i3p_i. Calculations 
for all these states are possible by means of our mesh HF method. However they are extremely tedious and time 
consuming and have not been accomplished in the present work. Indeed we will argue in the following that none of 
these states can be the ground state of the Li atom for intermediate field strength. 

It is quite evident that for the configurations containing a Is 2 pair of electrons the ls 2 3s lies higher in energy 
than the ls 2 2s configuration and that the ls 2 3p_i and ls 2 3d_2 configuration possess higher energy than the ls 2 2p_i 
configuration. Thus, the states with Is 2 pairs can be excluded from our argumentation of the ground state. Among 
the fully spin polarized configurations the levels of the configurations Is2p_i3s, Is2s3p_i, Is2s3<i_2, and Is2p_i3p_i 
are higher than that of the Is2s2p_i configuration (two components of the configurations are identical with those of 
Is2s2p_i and the third one is significantly higher). Thus from simple geometrical reasons only the Is2s3s configuration 
(mixed with the ls2s3do configuration) is a priori not excluded from becoming the intermediate ground state. In 
weak magnetic fields this state lies slightly lower than other doubly excited and autoionizing states and in this regime 
it is the lowest fully spin-polarized state. But the change of the ground state to the fully spin-polarized configuration 
takes place in the vicinity of 7 = 2 for which the 3s wave functions is much weaker bound than the 3d_2, 2p_i and 
even 2po orbitals. Due to this fact also the Is2s3s configuration can be excluded from becoming the ground state 
for any field strength. Indeed our calculations show that this state becomes higher in energy than the Is2s2p_i at 
7 w 0.16. 

Thus, the set ls 2 2p_i, Is J, 2s J, 2p_i |, and Is J. 2po J, 2p_i J, along with weak- ls 2 2s and strong-field Is j 2p_i J. 
3c?_2 4 ground states is comprehensive for the determination of the ground state of the Li atom in a magnetic field of 
arbitrary strength. 
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IV. NUMERICAL RESULTS 



The only work on the Li atom in a magnetic field with which we can compare our results is ref. In this 

reference HF calculations were performed for weak and intermediate magnetic field strengths. Table [j] contains the 
total energies obtained for the Li atom within our calculations in comparison with the data obtained in ]lT| . Our 
energy values coincide with those of ref. pd| ] for weak fields and lie substantially lower in the intermediate regime. At 
the upper boundary of the field region investigated in Jul the difference between and our energies is 0.0239 for 
the ls 2 2s state, 0.0205 for the ls 2 2p_i state, 0.0870 for the Is2s2p_i state, and 0.0458 for the ls2p 2p_i state. 

Our results on the total energies are illustrated in figures 1 and 2. These figures show in particular the ground 
state configurations for the different regimes of the field strength. One can conclude from table | and figures 1 and 
2 that the ls 2 2s configuration represents the ground state for < 7 < 0.17633, for 0.17633 < 7 < 2.153 the ground 
state configuration is ls 2 2p_i, and for 7 > 2.153 the ground state configuration is Is J, 2p_i J. 3<i_2 |. The state 
Is J, 2po I 2p_i I presented in |ll[] as the high field ground state appears not to be the ground state of the Li atom 
for any magnetic field strength. 

Figure 3 presents spatial distributions of the total electronic densities for the ground state configurations of the 
lithium atom. In each row these densities are presented for the limits of the corresponding field strength regions 
including the transition points and for some value of the intermediate field strength in between. For each separate 
configuration the effect of the increasing field strength consists in compressing the electronic distribution towards the 
z axis. For the Is2p_i3ci_2 configuration for which all single-electron binding energies increase unlimited for 7 — > 00 
a shrinking process of this distribution in z direction is also visible. For the ls 2 2p_i configuration this effect is not 
distinct for the relevant field strengths. For the ls 2 2s state the opposite effect can be observed: the 2s electronic 
charge distribution along the z axis expands slightly in weak magnetic fields. A characteristic feature of the transition 
points is an inflation of the electronic distribution in p direction during transitions from lower- to higher-field ground 
state configurations. This effect occurs due to the prevailing of the lowering in energy with changing quantum numbers 
(772 = to ui — — 1 for the transition point 7 = 0.17633 and S z = Ylu=i Sz m = — to S z = —3/2 for 7 = 2.153) over 
the raising of the energy due to more extended charge distributions in the p direction. 

The total binding energies of the configurations ls 2 2s, ls 2 2p_i, Is I 2s I 2p_i j, Is j 2p J, 2p_i j and Is J. 2p_i J. 
3d_2 I are presented in figure 4. These values do not include spin polarization terms and it can clearly be seen that 
the atomic ground state in a magnetic field does in general not possess the largest binding energy. 

Along with the total energy of the Li atom ground state we have obtained its ionization energies Ei dependent 
on 7. The total energy values of the ground state of the ion Li + are required for these calculations. The set of the 
ground state configurations of this two-electron ion is analogous to those of the helium atom |],[l(| and consists of 
the zero-field ground state Is 2 and the strong field fully spin-polarized state Is 1 2p_i !. Results of our calculations 
for these states are presented in table ||. The change of the ground state configuration takes place at 7 = 2.071814. 
Comparing tables [| and [n] one obtains the dependence of the ionization energy of the ground state of the Li atom on 
the magnetic field strength, as shown in figure 5. This curve exhibits three distinct points marked by dotted vertical 
lines. The first of them (from left to right) corresponds to the change of the ground state configuration of the lithium 
atom from ls 2 2s to ls 2 2p_i. The second corresponds to the change of the Li + ground state configuration from Is 2 to 
Is I 2p_i. And the third, very near to the second one, corresponds to the second change of the Li atom ground state 
configuration from ls 2 2p_i to Is J, 2p^i J. 3d_2 I- Table II provides the numerical data for the ionization energies. 
Tables | and || allow also obtaining ionization energies for other states presented in table |. 

In addition we show in figure 6 the total quadrupole moment 

Q zz = (#|3z 2 -r 2 |<&), r 2 =p 2 + z 2 (4) 

of different states of the atom as a function of the field strength. These dependencies illustrate the changes in the 
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density distribution of the electrons with increasing magnetic field strength. For weak and also to some extent for 
intermediate field strengths the main effect consists in compressing the wave function towards the z axis. This results 
in increasing Q zz values and a sign change of Q zz for the states with initially negative Q zz . For 7 > 10 the continuing 
compression towards the z axis practically does not affect Q zz due to the small values of (p 2 )- The values of Q zz 
decrease in this region for all the states considered with exception of the state Is J. 2p J, 2p_i j. This decrease of 
Q zz is associated with the decreasing value of (z 2 ) due to an increasing one-particle binding energy. For the states 
ls 2 2p_i and Is J, 2p_i I 3d_2 | all these binding energies become infinite for infinite strong fields. This results in 
Qzz —> as 7 — > 00. For the other states presented in Figure 6 at least one of the single-electron energies remains 
finite as 7 — > 00 and, in result, Q zz has a finite limit as 7 — > 00. 

V. SUMMARY AND CONCLUSIONS 

We have applied our 2D mesh Hartree-Fock method to a magnetized Li atom. The method is flexible enough 
to yield precise results for arbitrary field strengths and our calculations for the ground and several excited states 
are performed for magnetic field strengths ranging from zero up to 2.3505 • 10 8 T (7 = 1000). Our consideration 
was focused on the ground state of the Li atom. With increasing field strength this state undergoes two transitions 
involving three different electronic configurations. For weak fields up to 7 = 0.17633 the ground state arises from 
the field-free ls 2 2s configuration. For intermediate fields (0.17633 < 7 < 2.1530) the ground state is constituted by 
the ls 2 2p_i configuration and for 7 > 2.1530 the ground state configuration is the fully spin-polarized Is2p_i3c?_2 
configuration. We provide arguments which show that this configuration must correspond to the ground state in 
the strong field limit. Generalizing these arguments we could derive the high-field ground state configuration of 
arbitrary fully spin polarized atoms which are constituted by certain tightly bound hydrogen-like states. For example 
for atoms with six electrons (i.e. C and C-like ions) the high-field ground state is given by the fully spin polarized 
Is I 2p-i I 3cL 2 I 4/_ 3 I 5.g_4 J. 6/i_ 5 J. configuration. 

We have also calculated HF energies for the two Li+ ground state configurations Is 2 and Is J, 2p_i J.. The first 
of them forms the ground state at < 7 < 2.071814, the second one is the high-field ground state configuration for 
7 > 2.071814. These calculations allowed us to obtain the Li atom ground state ionization energy Ej dependent on 
the magnetic field strength. This dependence, opposite to the analogous dependence for the total and binding energies 
is not monotoncous and contains both areas of increasing values of Ej and a domain of decreasing behaviour between 
7 = 2.071814 and 7 = 2.1530. Furthermore we have studied the quadrupolc moment of the atom and show how its 
complicated behaviour with changing field strength can be explained through the field dependence of the different HF 
orbitals. 

Two remarks are in order. Our HF results do not include the effects of correlation. To take into account the latter 
would require a multiconfigurational approach which goes beyond the scope of the present paper. We, however, do 
not expect that the correlation energy changes our main conclusions like, for example, the transitions in the ground 
states configurations or the behaviour of the ionization energies depending on the field strength. With increasing field 
strength the effective one particle picture should be an increasingly better description of the wave function and the 
percentage of the correlation energy should therefore decrease. For the case of hydrogen it is well-known that in the 
high field regime (7 >> 10 2 ) mass correction terms due to the finite nuclear mass become relevant i.e. are no more 
negligible in comparison with the Coulomb binding energies. The most important mass corrections can be included 
by replacing the electron mass through its reduced mass and results from the infinite nuclear mass calculations are 
related to those with the reduced mass via a scaling relation. In the case of the much heavier Li atom these effects 
are expected to be much smaller. 

Apart from the Li atom other species i.e. three-electron objects are expected to be in particular of astrophysical 
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interest: the three-electron ions formed by the nuclei He, C, O, and Nc possess a high abundance in the universe. To 
study these systems is the subject of a separate investigation. 
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Figure Captions 

Figure 1. Total energies of the Li atom as a function of the magnetic field strength (solid lines marked by centered 
symbols). Dotted lines are energies of two electronic configurations of the Li + ion: (a) low- field ground state Is 2 ; (b) 
high- field ground state ls2p_i. 

Figure 2. The same as in figure 1 in the relevant regime of transitions of the ground state configurations. 

Figure 3. Contour plots of the total electronic densities for the ground state of the Li atom. The densities for 
neighbouring lines are different by a factor of e. 

Figure 4. Binding energies of various states of the Li atom as a function of the magnetic field strength. 

Figure 5. Li atom ground state ionization energy Ej for a broad range of field strengths. Transition points are 
marked by broken vertical lines. The first transition (from left to right) corresponds to the change of the ground 
state configuration from Is 2 2s to ls 2 2p_i. The second transition corresponds to the change of the Li + ground state 
configuration from Is 2 to ls2p_i. Third transition of the Li ground state configuration from ls 2 2p_i to Is2p_i3d_2. 

Figure 6. Quadrupole moment of the Li atom depending on the magnetic field strength. 
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TABLE I. Total energies of several electronic ground and excited states of the Li atom in the regime of field strength 
7 = 0, ...,1000 
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0.0020 


-7.43375 






-7.36709 








-5.36288 








-5.23586 






-5.08976 


0.0050 


-7.43522 






-7.37002 








-5.36884 








-5.24182 






-5.09852 


0.0090 


-7.43713 


-7.4371 


-7.37387 


-7.3738 


-5.37673 


-5.3767 


-5.24973 


-5.2497 


-5.10988 


0.0100 


-7.43760 






-7.37481 








-5.37871 








-5.25170 






-5.11268 


0.0180 


-7.44125 


-7.4412 


-7.38218 


-7.3832 


-5.39429 


-5.3943 


-5.26734 


-5.2673 


-5.13433 


0.0200 


-7.44214 






-7.38397 








-5.39817 








-5.27121 






-5.13960 


0.0500 


-7.45398 






-7.40844 








-5.45442 








-5.32786 






-5.21281 


0.0540 


-7.45537 


-7.4553 


-7.41141 


-7.4114 


-5.46168 


-5.4617 


-5.33521 


-5.3352 


-5.22199 


0.1000 


-7.46857 






-7.44176 








-5.54149 








-5.41643 






-5.32140 


0.1260 


-7.47408 


-7.4739 


-7.45650 


-7.4565 


-5.58376 


-5.5837 


-5.45992 


-5.4599 


-5.37371 


0.17633 


-7.48162 






-7.48162 
























0.1800 


-7.48204 


-7.4814 


-7.48330 


-7.4832 


-5.66585 


-5.6656 


-5.54555 


-5.5455 


-5.47568 


0.2000 


-7.48400 






-7.49220 








-5.69451 








-5.57585 






-5.51151 


0.5000 


-7.47741 






-7.58790 








-6.04787 








-5.96957 






-5.97052 


0.5400 


-7.47351 


-7.4731 


-7.59709 


-7.5965 


-6.08746 


-6.0844 


-6.01603 


-6.0159 


-6.02414 


0.9000 


-7.42504 


-7.4240 


-7.65628 


-7.6563 


-6.40175 


-6.3993 


-6.39613 


-6.3956 


-6.46061 


1.0000 


-7.40879 






-7.66653 








-6.48029 








-6.49248 






-6.57081 


1.2600 


-7.36226 


-7.3609 


-7.68288 


-7.6820 


-6.67494 


-6.6720 


-6.72931 


-6.7284 


-6.84122 


1.8000 


-7.24603 


-7.2446 


-7.67657 


-7.6747 


-7.05430 


-7.0403 


-7.17326 


-7.1711 


-7.34723 


2.0000 


-7.19621 






-7.66246 








-7.18889 








-7.32494 






-7.52003 


2.071814 


-7.17745 






-7.65600 








-7.23650 








-7.37799 






-7.58047 


2.1530 








-7.64785 






















-7.64785 


2.1600 








-7.64711 


-7.6459 


-7.29445 


-7.2826 


-7.44218 


-7.4404 


-7.65361 


2.5000 


-7.05619 






-7.60351 








-7.51255 








-7.71826 






-7.92532 


3.0000 


-6.89559 






-7.51516 








-7.81834 








-8.00837 






-8.29920 


3.6000 


-6.67874 


-6.6640 


-7.37638 


-7.3627 


-8.16336 


-8.1159 


-8.37214 


-8.3564 


-8.71464 


3.9600 








-7.27826 


-7.2722 


-8.35994 


-8.3165 


-8.57739 


-8.5578 


-8.94929 


4.3200 








-7.17026 


-7.1655 


-8.54941 


-8.5075 


-8.77415 


-8.7526 


-9.17442 


4.6800 








-7.05326 


-7.0391 


-8.73233 


-8.6767 


-8.96327 


-8.9371 


-9.39099 


5.0000 


-6.08811 






-6.94230 








-8.88981 








-9.12554 






-9 57694 


5.0400 








-6.92800 


-6.9050 


-8.90918 


-8.8375 


-9.14546 


-9.1160 


-9.59977 


5.4000 


-5.90113 


-5.8772 


-6.79517 


-6.7747 


-9.08045 


-9.0035 


-9.32134 


-9.2755 


-9.80147 


7.0000 


-5.08909 






-6.12670 








-9.78357 








-10.03896 






-10.62578 


10. 


-3.35777 






-4.61777 








-10.91059 








-11.17886 






-11.93902 


20. 


3.49120 






1.70565 








-13.69420 








-13.96582 






-15.16260 


50. 


27.6916 






24.97942 








-18.8012 








-19.0436 






-21.0505 


100. 


71.807 






68.1735 








-23.987 








-24.1946 






-27.0192 


200. 


164.371 






159.5749 








-30.559 








-30.7327 






-34.5850 


500. 


451.69 






444.9033 








-41.821 








-41.959 






-47.5583 
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1000. 



939.54 



930.84308 



52.65 



-52.771 



-60.0589 
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TABLE II. Energies of the low- and high-field ground states of the ion Li + and the ionization energy of the ground state of 
the Li atom Ei for field strengths 7 = 0, 1000. 



7 


Is 2 


ls2p_i 


E,{U) 


0.0000 


-7.23642 


-5.02469 


0.19633 


0.0010 


-7.23642 


-5.02619 


0.19684 


0.0020 


-7.23642 


-5.02769 


0.19733 


0.0050 


-7.23641 


-5.03218 


0.19881 


0.0100 


-7.23641 


-5.03963 


0.20119 


0.0200 


-7.23639 


-5.05442 


0.20575 


0.0500 


-7.23623 


-5.09797 


0.21775 


0.1000 


-7.23567 


-5.16789 


0.23290 


0.17633 


-7.23411 


-5.26874 


0.24751 


0.2000 


-7.23345 


-5.29873 


0.25875 


0.5000 


-7.21798 


-5.64006 


0.36992 


1.0000 


-7.16401 


-6.11462 


0.50252 


2.0000 


-6.96300 


-6.89408 


0.69946 


2.071814 


-6.94440 


-6.94440 


0.71160 


2.1530 


-6.92278 


-7.00057 


0.64729 


2.5000 


-6.82347 


-7.23258 


0.69275 


3.0000 


-6.66237 


-7.54672 


0.75248 


5.0000 


-5.85051 


-8.62943 


0.94751 


7.0000 


-4.84725 


-9.52492 


1.10086 


10. 


-3.11092 


-10.65131 


1.28771 


20. 


3.74896 


-13.42974 


1.73286 


50. 


27.96465 


-18.52548 


2.5250 


100. 


72.09337 


-23.69994 


3.3193 


200. 


164.66867 


-30.26077 


4.3242 


500. 


452.0032 


-41.50393 


6.0544 


1000. 


939.87976 


-52.3230 


7.7359 
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